Unwinding relaxation dynamics of polymers 
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The relaxation dynamics of a polymer wound around a fixed obstacle constitutes a fundamental 
instance of polymer with twist and torque and it is of relevance also for DNA denaturation dy- 
namics. We investigate it by simulations and Langevin equation analysis. The latter predicts a 
relaxation time scaling as a power of the polymer length times a logarithmic correction related to 
the equilibrium fluctuations of the winding angle. The numerical data support this result and show 
that at short times the winding angle decreases as a power-law. This is also in agreement with 
the Langevin equation provided a winding-dependent friction is used, suggesting that such reduced 
description of the system captures the basic features of the problem. 

PACS numbers: 82.35.Lr, 36.20.Ey, 61.25.hp 
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The dynamics of polymers subject to spatial or topo- 
logical constraints has received quite some attention in 
recent years. Interesting examples are the translocation 
of DNA from a narrow pore (for a recent discussion see 
e.g. [l| and references therein) or the dynamics of su- 
percoiled DNA (see e.g. Q). An important question is 
whether the complex polymer dynamics can be described 
by a simple equation of motion, using a one-dimensional 
reaction coordinate. This issue arises, for instance, in the 
context of polymer translocation (from a pore in a wall) 
where it was shown that the Langevin equation fails to 
reproduce simulation results Q. This failure motivated 
extensive studies. Various models were put forward, as 
the generalized Langevin equation with a memory ker- 
nel jj] , or a deterministic two-phase model [B[ . 

The aim of this Letter is to study analytically and nu- 
merically the unwinding relaxation dynamics, which also 
belongs to the above class of problems. The equilibrium 
winding angles for polymers were intensively studied in 
the past [6|48|. These studies are relevant for a series of 
problems in physics, as e.g. models for the behavior of 
flux lines in high-Tc superconductors Q . The relaxation 
dynamics of unwinding has been much less studied jilfioj]. 
though it is a problem of relevance in DNA melting dy- 
namics, but also as a fundamental issue of polymer dy- 
namics involving twist and torsion. 

We consider a polymer initially wound around a long 
impenetrable bar (see Fig. [T]), to which it is attached at 
one end. Since this is an entropically highly unfavorable 
situation, the polymer will unwind, starting at the loose 
end; and given enough time, it will relax towards the 
equilibrium state in which it is no longer winding around 
the bar. To monitor the unwinding process, we keep track 
of the winding angle of the last monomer of the poly- 
mer, which measures the angle accumulated by the chain 
around the bar from the first attached monomer to the 




FIG. 1. (Color online) Snapshot of a polymer configuration on 
the fee lattice, during the unwinding from a bar. The winding 
angle is defined as = ^ =1 A^i, where Aflj^+i is the 
difference in angles between monomer i + 1 and i measured 
with respect to the bar. The hue follows the monomers order 
from i = l ("O", attached to the bar, blue online) to i — L 
( "E" , red online) . The configuration displays a tightly bound 
helix (OB), a loose helix (BA), and a free end (AE). 



last free one. We treat the case of polymers with internal 
excluded volume by studying a self-avoiding walk (SAW) 
and we support our arguments by also investigating the 
motion of a random walk (RW). Compared to the more 
complex unwinding of a double stranded DNA helix, the 
advantage of dealing with a single polymer around a fixed 
obstacle is that this winding angle provides a well-defined 
"reaction coordinate" . 

The numerical calculations were performed using lat- 
tice polymers, specifically L-step RWs on a square lat- 
tice and SAWs on a face-centered-cubic (fee) lattice. 
An update consists of a local corner flip or an end-flip 
move (Rouse dynamics), and a time step includes L up- 
dates at random locations. The initial configuration for 
a RW is constructed by the repetition of a sub-walk 
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FIG. 2. (Color online) Simulations (solid lines) of average 
winding angles vs. time in a semi-logarithmic plot for (a) RW 
in two dimensions and (b) SAW in three dimensions. Only 
data for two polymer lengths are shown (see legends). Dashed 
lines are fit of the exponential decay at long times. Inset of 
(a): zoom-in of the short-time behavior with extrapolated 
intercepts 6i n . Inset of (b): plot of 9i n vs. polymer length L, 
for the SAW data, showing a linear behavior ~ L. 



of 8 monomers winding around the (0,0) site (repre- 
senting the bar), with a resulting initial winding angle 
#o = tiL/A. Similarly, for the SAW on the fee lattice we 
repeat a helix formed by 6 steps {8q = irL/3) around the 
bar in the direction (1, 1, 0). Figure[3]shows a plot of 8 vs. 
time in a semi-logarithmic scale obtained from numerical 
simulations: one distinguishes a long-time regime where 
8 relaxes exponentially and a short-time regime that de- 
viates from the exponential decay. We will discuss the 
two cases separately. 

Our analytical scheme is based on a one-dimensional 
Langcvin equation for the variable 8: 



dd dT{8,L) 
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with J- the equilibrium free energy for a polymer of length 
L and winding angle 8, j T the torque friction and rj a 
noise term. The main focus is on the time evolution 
of the average winding angle (8) (indicated with 8 for 
simplicity) and not in fluctuations, so the noise term will 
be neglected. 

The degrees of freedom parallel to the bar are not rel- 
evant for a RW and one can restrict the study to a two- 
dimensional walk, with the bar replaced by an excluded 
site. For a planar RW the free energy is known exactly @ : 



T RW (8,L) = -k B T\og 



cosh 



logL 
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For a SAW wound around the bar recent numerical sim- 
ulations suggest a similar scaling form [ll| : 



T S aw{0,L) = -k B T log 



V 



(logL)O- 



(3) 



where p() is the probability distribution of winding angles 
obtained from equilibrium Monte Carlo sampling. Here 



the exponent 0.75 is a numerical estimate Since 
both free energies involve a scaling variable 8/(\ogL) a , 
with a = 1 for RWs and a sw 0.75 for SAWs, we can 
analyze the two processes on equal footing. 

We focus first on the longest relaxation time. Eqs. © 
and ^ are quadratic for small 8. Hence using the lowest- 
order term and neglecting other proportionality factors, 
for small angles one obtains an equation of the form 
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dt 7 T (logL) 
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In order to gain some insights on the L-dcpcndence of the 
torque friction 7 T one can consider a particle rotating 
at a fixed distance R from an origin and subject to a 
constant tangential force /. The Langevin equation in 6 
is of the form 7 T 4f = = fR, where f2 is the torque. The 
equation can be transformed into a cartesian coordinate 
x = R8, yielding 



/ = 



7 r Rd8 7 T dx 
~R2~dt 



dx 

= l¥~dt = 1 ~dl ' 



(5) 



where 7 is the friction associated with linear displace- 
ment. This implies that 7 T = R 2 ^. By integrating 
over L monomers we obtain 7 ~ L and hence an av- 
erage torque friction 7 T ~ L 1+2v , where the Flory expo- 
nent v describes the average end-to-end squared distance 
(R 2 ) ~ L 2v for a polymer in equilibrium: v — 1/2 for a 
RW, while v ~ 0.588 for a three-dimensional SAW. Plug- 
ging the estimated 7 r into Eq. ^ one finds the following 
relaxation time-scale 



tl 



- l+2u 



(logL) 
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If hydrodynamic effects arc included, the friction for 
linear displacement grows as 7 ~ V ', and the relax- 
ation time becomes tx ~ L 3u (logL) 2a . Note that the 
leading term of Eq. (J6j> is similar to the Rouse time 
T Rouse ^ i}+ 2v t which is the equilibration time of a free 
polymer 12j . This is also a lower bound for the unwind- 
ing relaxation time, i.e. 77, > r|^ ousc , as the attachment 
to the bar and its stcric hindrance arc unlikely to speed 
up the equilibration process. 

In the simulations we determined the total unwinding 
time r£, i.e. the average time needed for the unwinding 
process to be completed. We defined it as the time it 
takes to reach 8 = for the first time. As the polymers 
are initially wound to 8q ~ L, one has to take into ac- 
count that the relaxation starts from a higher winding 
angle for longer polymers. The analysis of the numer- 
ical data (see Fig. [2]) shows that the asymptotic decay 
is well-fitted by 8{t) = 8i n exp(— t / 'tl) and the intercept 
6 i n scales linearly with L. Hence the condition 8 (r£) ~ 1 
gives 



r L logL (logL) 



\2a + l 



(7) 



Thus, the total unwinding time r£ differs by a factor 
log L from the relaxation time-scale tl ■ 



3 



20 



* 15 

2 10 



1 1 

(a) 

_ 






- 2.65(2)' 




- 2.32(2) 


r , 1 


i 


2 4 


6 8 



logL 




id 

CO 



2.5 





(c) 

/ 


1 1 1 1 1 \y 1 

(d) /j^ 


/ — S(aL),a=6.5 

/ i \ 1 1 1 1 'i 


§(aL),a=4.5 



0.1 0.2 0.3 0.1 0.2 0.3 0.4 
1 / log L 1 / log L 

FIG. 3. (Color online) (a) Plot of logr^ vs. logL for the RW, 
with lengths up to L = 2048. Averages are from 10 4 inde- 
pendent runs for L < 200 down to 200 for L = 1536, 2048. 
Dashed and solid lines are fits in the short and long L regimes, 
showing a systematic variation in the exponent, (b) Same as 
(a) for SAWs, with lengths up to L = 144 (10 s independent 
runs per L). (c) Squares: Plot of the RW running expo- 
nent 8(L) obtained from simulations [estimated by a centered 
difference of data in (a)] as a function of 1/logL. Dashed 
and solid lines represent the scaling 5(L) and S(auwL), with 
o-rw = 6.5, from Eq. (5), respectively, (d) The same for 
SAWs data from panel (b), with asAw = 4.5. 



Plots of logr£ vs. logL are shown in Fig. [3ja) (RW) 
and[3Jb) (SAW). In order to analyze the data appropri- 
ately we computed 8(L), defined as the "local" slope in 
the log t£ vs. log L plot for a given size L. Eq. (JT]) implies 



5(L) 



dlogL 



= l + 2v 



2a + 1 
logL 



(8) 



Figures 03c) and (d) (squares) show the numerical es- 
timates of 8(L) vs. 1/logL for the RW and SAW, re- 
spectively. The asymptotic scaling predicted by Eq. J5]) 
implies a straight line for 8(L) when plotted as a func- 
tion of 1/logL (dashed lines in Figs. |3Jc) and (d)). 
The curvature in the data indicates that further finite- 
size corrections should be included. In order to ra- 
tionalize them we introduce a finite-size scaling ansatz 
8{aL) = 1 + 2v + (2a + 1)/ log(aL) in which an ampli- 
tude "a" is included in the logarithmic factor as a single 
fitting parameter. The best fit of 8(aL) to the data points 
produces the two solid lines in Figs. [3^c) and (d). The 
data for the RW (c) are in excellent agreement with the 
ansatz, while the SAW data are less conclusive: they in- 
volve much heavier computations and are thus restricted 
to much shorter polymers. 

We consider next the early-time dynamics of 9. Fig- 
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FIG. 4. (Color online) 6 -8(t) versus t for the RW (a) and the 
SAW (b). The log-log scale highlights the short-time regime, 
which behaves as a power law (Eq. ©) with p ~ 0.44 (RW) 
and p~ 0.43 (SAW). 



ureS] shows a plot of 9 — 9(t) vs. t in log-log scale. For 
t < 10 4 the data are fitted by a power-law behavior 



9(t) ~ t" . 



(9) 



with p w 0.44 and p w 0.43 for RW and SAW, respec- 
tively. To understand this behavior we consider again 
the Langevin equation ((T|). At very high winding, where 
6(t) ~ 8q ~ L the torque due to free energies in Eqs. ([2]) 
and 43} arc of little use as they describe equilibrium fluc- 
tuations for small 8's. In the early stages of the dynamics 
we expect unwinding only near the free end, regardless 
of the polymer length. The decrease of 9 is then linearly 
related to the length of the unwound part of the poly- 
mer, and to leading order also linearly to the increase 
in entropy. We thus assume that the torque is constant 
(L- independent), tq = —^jf = const. At high winding 
the friction decreases, as the part of the polymer which 
is tightly wound around the bar does not contribute to 
it. For 9q — 9 <C L, the friction coefficient should depend 
only on the difference 9q — 9. Let us consider a friction 
coefficient vanishing as a power-law as j T (9) ~ (#o — ®Y ■ 
Combining this ansatz for j T with the argument for a 
constant torque tq, from ([T]) one obtains 



m\ dd 



9 -9r-(6 -9)~T 



(10) 



which integrated in time, and using the initial condition 
9(0) = 8q, yields a power-law scaling as that given in 
Eq. (0 with p= 1/(1 + a;). 

To estimate the exponent x we introduce two differ- 
ent types of hypotheses about the shape of the polymer 
in the early stages of unwinding. These are sketched in 
Fig. [5ja) and (b). In the case (a), we consider a tightly 
wound polymer for a length L — I and an unwound loose 
part of length I and assume that the latter is equili- 
brated. We denote the winding per unit length in the 
wound part with Aa>i and that of the loose part with 
Aw 2 (Aw 2 ~ in the case of Fig.[5[a)). The winding an- 
gle 9 = (L — Z)Awi = O — iAwi, hence 9 -9 = lAui. As 
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FIG. 5. Two possible configurations of polymers during un- 
winding, (a) A tight helix of length L — I connected to a loose 
equilibrated end of length I. (b) The part of the polymer 
detached from the bar here has still some winding. The expo- 
nents p governing the early-time decay of the winding angle 
predicted in the two cases are given. 



shown above in the discussion of the late-time relaxation, 
an equilibrated polymer of length I has a torque friction 
scaling as l 1+2v , therefore 7 T (0) ~ (0o — 0) l+2y which im- 
plies p = l/(2v + 2) = 1/3 for a RW and p = 0.31 for a 
SAW. An alternative conformation is shown in Fig.[5jb). 
In this case we consider a "looser" helix of length I with 
density of winding per unit length Acl>2 > connected to 
a tightly wound helix of length L — I. Only the former 
contributes to the friction. In addition we assume that 
the looser helix does not change its radius and pitch in 
time (thus Au>2 is constant). This seems reasonable at 
least for the early times of the dynamics. We then have 
= (L - Z)Awi + IAuj 2 = 0o - l(Au>i - Aw 2 )- As the 
loose helix maintains its shape while growing the friction 
is simply proportional to its length: r y T {9) ~ I ~ (0o — 0), 
which yields p=\j1 both for a RW and a SAW. 

The conformations of Fig. [5] are of course "idealized" 
and should represent two extreme cases. In (a) the loose 
end stretches out from the bar causing a more rapid in- 
crease in the friction compared to (b). The exponent 
p = 1/2 predicted for the case (b) is quite close to 
p « 0.45 found in simulations. Snapshots such as that in 
Fig. [T] suggest that the actual polymer conformations are 
hybrids of those in Fig. [5J Starting from the free end, one 
notices a very loose part which does not add much to the 
winding angle (segment AE in Fig[T]). This is reminis- 
cent of the loose equilibrated end of Fig. [Ha). There is 
then an intermediate part (BA in Fig. [1]) wound around 
the bar, but not tightly, resembling the loose helix of 
Fig-ETb). 

In conclusion, in this Letter we investigated numeri- 
cally the relaxation dynamics of polymers wound around 
a fixed obstacle and we have provided an analytical 
scheme based on a Langevin equation for the winding an- 
gle. Studying such equation in the late relaxation stage 
we predict the scaling form of the friction and conse- 
quently of the unwinding time-scale, which involves log- 
arithmic corrections to the power-law of the chain length. 



The same equation is also useful in the regime at short 
times, where a friction depending on the unwinding is 
needed to describe the observed scaling of the winding- 
angle . The two cases analyzed numerically, a SAW and 
a RW, provide a consistent picture of the dynamical be- 
havior. Although logarithmic factors are notoriously dif- 
ficult to study in simulations, finite-size scaling extrap- 
olations of our results are compatible with the predic- 
tions of the Langevin equation. It is possible that such 
strong corrections affect also the unwinding of two poly- 
mers from a double-helical conformation. A recent nu- 
merical study (lo| yields an unwinding time scaling as 
t* l ~ L 2 - 58 ; numerically, this scaling is consistent with 
that of the running exponent found in this work for the 
longest polymers (see Fig. 05b)). It is thus plausible that 
the relaxation time of an unwinding double-helix is also 
captured by Eq. ((7j). Besides delving new fundamental 
aspects of polymer dynamics and providing a reference 
case for DNA denaturation dynamics, this study may also 
serve as a basis for other types of investigations involv- 
ing rotational dynamics, as for instance the relaxation of 
plectonemic structures which form in overtwisted DNA. 
Modeling the statics and dynamics of DNA plcctonemcs 
has been of recent great interest [H, [l3l - [l5l | . 
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